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We present a detailed theory of spectacular semiclassical catastrophes happening during the 
time evolution of a kicked quantum rotor (Phys.Rev. Lett. 87, 163601 (2001)). Both two- and 
three-dimensional rotational systems are analyzed. It is shown that the wave function of the rotor 
^/-^ | develops a cusp at a certain delay after a kick, which results in a sharply focused rotational wave 

f"^ ■ packet. The cusp is followed by a fold- type catastrophe manifested in the rainbow- like moving 

' angular singularities. In the three-dimensional case, the rainbows are accompanied by additional 

CNj , singular features similar to glory structures known in wave optics. These catastrophes in the time- 

<— ( . dependent angular wave function are well described by the appropriate tools of the quasiclassical 

wave mechanics, i.e. by Airy and Bessel approximations and Pearcey's functions. A scenario of 
K— j , "accumulative squeezing" is also presented in which a specially designed train of short kicks produces 

■ an unlimited narrowing of the rotor angular distribution. This scenario is relevant for the molecular 

| alignment by short laser pulses, and also for atom lithography schemes in which cold atoms are 

focused by an optical standing wave. 
7-H ' PACS numbers: 32.80. Lg, 03.65. Sq, 05.45.-a 

> : 
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I. INTRODUCTION 

\o : 

The driven rotor is a standard model for classical and quantum-mechanical nonlinear studies. The interest to the 
CO ■ subject has substantially increased because of new laser experiments on molecular alignment (orientation) (for a 
recent review see [1]), and due to the quantum optics realization of the 5-kicked rotor [2,3]. Alignment of molecules 
by using laser pulses that are much shorter than the typical molecular rotational times is of special interest. At 
O ,' these conditions, a coherent rotational wave packet is generated that takes an angularly squeezed shape after the 
pulse is over [4-15]. This phenomenon opens new prospects for numerous applications in chemical physics and non- 
linear optics which require enhanced molecular orientation (alignment) for a limited time but at field-free conditions. 
Generation of ultra-short laser pulses [16,17] and control of high harmonics generation [18] are only few examples to 
mention. 

Motivated by these reasons two of us have revisited the problem of 5-kicked quantum rotor recently [19]. It was 
shown that free evolution of a strongly kicked rotor is subject to several semiclassical catastrophes with course of 
time, which are manifested in the angular distribution singularities. In particular, a cusp develops in the angular wave 
function at certain time after the pulse, providing an extremely narrow peak in the direction of the kick. This effect 
is a time-domain analog of the focusing phenomenon in wave optics. Moreover, the subsequent spread of the focused 
rotational wave packet demonstrates a wave front fold leading to the moving rainbow-type singularities in the angular 
distribution. These effects can be observed both in two-dimensional and three-dimensional rotational systems. The 
2D rotational model is important not only due to its relative simplicity, but also because it describes the interaction 
of an ensemble of cold atoms with sinusoidal standing light waves. In 3D case the geometry of the wave bifurcations 
is even richer [20]. In particular, the angular cusp (focus) transforms not only to the moving rainbow feature, but 
also to a persistent singular peak in the direction of the kick. A little later such a peak appears also in the opposite 
direction. These two peaks are reminiscent of the forward and backward glories well known in the scattering of light 
by small droplets [21]. As time goes on, the rotational wave packet exhibits a rich scenario of fractional revivals [22] 
with the presence of many angular rainbows and foci. 

The cusp, fold and glory are well known catastrophes in wave optics and have been studied since long ago (see, 
e.g. [23,24] and references therein). A complete study of the morphologies of caustics in optics has been given in 
[25]. Catastrophes are also known in the field of atomic and nuclear physics where some of them (like rainbow and 
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glory) were studied at a very early stage [26]. Besides many experimental works to unravel caustics [27,28], a large 
number of theoretical studies were devoted to the semiclassical description of the occurring singularities [29-31,25,32]. 
The simplest semiclassical approaches are based on the stationary phase approximation for integral representations 
of the related wave functions. They can be successfully used if the stationary phase points are well separated from 
each other. However, the semiclassical catastrophes like rainbow and glory involve different kinds of non-isolated and 
merging stationary points. Various uniform approximations have been developed for such a case [33-36]. 

In the present study we apply semiclassical methods to the time-domain cusp, fold and glory phenomena exhibited 
by rotors subject to strong kicks. At the beginning of Section II we specify two different kinds of coupling to the 
external field (dipole-type and polarization-type interaction). For pedagogical reasons, we then continue with the 
classical analysis of the angular distribution of an ensemble of kicked rotors. We start with the 2D case, and consider 
additional effects specific to three-dimensional geometry after that. The same approach is kept in the following 
parts of Section II, where a fully quantum description of the problem is given. Section III presents various integral 
representations for the wave function of a kicked quantum rotational system. The next three Sections are devoted 
to semiclassical theory of the time-domain focusing, rainbow and glory effects for the kicked rotational systems, 
respectively . In Section VII we provide a classical analysis of the time-dependent kicked rotor angular distribution at 
finite temperature. Section VIII considers a systematic squeezing of the rotor angular distribution by a train of multiple 
pulses. In [19] we proposed a multi-pulse excitation scheme that exhibits an accumulative squeezing when increasing 
the number of kicks. After that, different strategies leading to enhanced squeezing have been suggested, compared 
and optimized for an ensemble of cold atoms in a pulsed optical lattice [37] and molecular systems [15]. Recently, 
accumulative squeezing was demonstrated experimentally in the atom optics system [38] . Section VIII provides details 
of the asymptotic behavior of the accumulative squeezing scenario in the quasiclassical regime. Finally, the results of 
the paper are summarized in Section IX. The paper has also two Appendices with details of calculations. 



II. KICKED ROTOR SYSTEM 



Under certain conditions, the process of molecular orientation (or alignment) by laser fields can be described by a 
strongly driven 3D rigid rotor model. The Hamiltonian of a 3D driven rotor is 

H = ^ + v(e,t), (i) 

where L is the angular momentum of the rotor, and / is its moment of inertia. For a linear molecule, having a 
permanent dipole moment [i, and driven by a linearly polarized field, the interaction potential is 

V(0,t) = -fj£(t) cos(0) (2) 

where £{t) is the field amplitude (i.e., of a half-cycle pulse) and 9 is the polar angle between the molecular axis and 
the field direction. In the absence of interaction with a permanent dipole moment, the external field couples with the 
induced molecular polarization. For nonresonant laser fields, this interaction, averaged over fast optical oscillations, 
is (see, e.g., [39,40]) 

V{6 1 t) = -- A £ 2 {t)[(a\\ -a ± )cos 2 (d)+a ± }. (3) 

Here ay and a± are the components of the polarizability, parallel and perpendicular to the molecular axis, and £(t) 
is the envelope of the laser pulse. The quantized motion of a rigid rotor is convenient to describe using dimensionless 
time scale, t = th/I. Introducing also the dimensionless interaction strength e p = ^i£(t)I/H 2 and ej = l/4£ 2 (i)[(a|| — 
a±)I/h 2 , the Hamiltonian Eq.(l) can be written either as 

H = ^-£ p (t)cos0, (4) 

or as 

L 2 

H= — -6 1 {t)cos 2 6, (5) 

for the two above interaction types. 
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In this work, we consider laser pulses that are so short that the motion of the rotor is " frozen" during the interaction 
with the pulse (S- kicked rotor). In this case, the action of the laser field is characterized by the time- integrated 
interaction intensity (kick strength) defined as 




Although we started our discussion with a three-dimensional rotor, it is equally interesting to study also the related 
two-dimensional systems. This is justified both pedagogically and by the existence of physical systems (i.e. cold 
atoms in standing optical waves [37,38]) that are described by such a model. 



A. Classical description of an ensemble of 2D kicked rotors. 

We start with the classical dynamics of an ensemble of free 2D rotors having a permanent dipolc moment (see 
Eq.(4)), which are kicked by a short pulse. The second interaction type (Eq.(5)) will be treated briefly after that. 

Consider a collection of motionless rotors subject to a pulse applied at r = 0. A particular rotor located initially 
at an angle 9o will be found at the same angle just after the kick, but having the angular velocity —Psm(9o). The 
time-dependent position of this rotor at time r is given by 

9 = O - Prsin((9 ) (mod27r). (6) 

For rotors starting their motion in the region of small angles (9 « 1), the acquired velocity is linearly proportional 
to the initial angle. As the result, all such rotors will arrive at 9 = ("focal point") at the same time Tf. 

Tf = VP 

providing a sharp peak in the angular distribution. This phenomenon is quite similar to the focusing of light by a 
thin optical lens. Certainly, the angular focusing is not perfect because of the aberration mechanism (deviation of the 
cos(9) potential from the parabolic one). For r < r/ , Eq.(6) presents a one-to-one mapping 9(9q) (see Fig. 1 (a)). At 
r = Tf the curve 9(9o) touches the horizontal axis (see Fig. 1 (b)), i.e. the mapping turns degenerate. At r > 17 the 
angle 9 becomes a multi- valued function of 9 (Figs. 1 (c) and 1 (d)). The probability of finding a rotor in a certain 
angle element d9 is determined by the initial distribution function, /o(#o) as follows: 

f(6)d6 = Y,f (6Z)\d6Z\. (7) 

a 

The summation in Eq.(7) is performed over all branches of the function 6 0(6) defined by (6). The modulus sign is 
needed to cover the branches in which 9 is a decreasing function of 9o- The time-dependent angular distribution of 
the ensemble of rotors is then given by 

It follows immediately from Eq.(8) that even for a smooth initial distribution, the function f(9,r) exhibits a singular 
behavior near the angles where d9/d9 — 0. This condition (together with Eq.(6)) defines the critical initial angle 
cos 9q = 1/(Pt) for which two classical trajectories coalesce at the same final angle 9 r (see Fig. 3 (a)). We may call 
such an event a rainbow-like catastrophe in analogy with the similar optical phenomenon. As is clear from Fig. 2, 
two such rainbows appear at r = r/, and they persist forever after that. Fig. 2 shows focus and rainbow singularities 
for the uniform initial angular distribution. The time-dependent rainbow angle 9 r is given by 



9 r = - arccos(l/PT) + ^(Pt) 2 - 1 

The other rainbow can be found at — 9 r . 

In the case of polarization- type interaction Eq.(5) we simply replace Eq.(6) by 

9 = 6» - Pt sin(20 o ) (mod 2tt) 

Similarly to the previous case, the focusing happens at Tf = 1/(2P). However, it can be easily seen that the focusing 
occurs simultaneously at 9 = and 9 — ir now. Because of the symmetry of the angular potential, two pairs of 
rainbows appear in the angular distribution at r > 17 . 
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B. Classical description of an ensemble of 3D kicked rotors. 



We use spherical coordinates in which field direction is taken as z-axis, and the angles 9, ip are limited by ^ 9 ^ n, 
^ ip $J 2ir. We restrict ourselves by axially symmetric distributions only. For the dipole-type coupling, the 
trajectories of the kicked rotors are obtained from Eq.(6) using the above restrictions for 9 and 9o (see Fig. 3 (b)). 
The time-dependent distribution function obeys the rule similar to (8): 

f(e T)= y /(flg,r = 0)sm(flg) 

Here the sum runs again over the branches of the function 9q(9) shown at Fig. 3. Clearly, we obtain focusing at r = r/ 
and a rainbow- like singularity at t > r/ (whenever \d9/d9 \ a = 0). The rainbow first appears at the north pole (9 = 0) 
and takes a shape of a ring moving on the sphere toward the south pole. The rainbow ring shrinks to a point at time 
t'j defined from 



arcC os(pL) + ^{p T ' f Y 



After that, the rainbow changes the direction of motion and runs toward the north pole. These cycles repeat over 
and over again. 

Comparing Eq.(8) and Eq.(9), we may identify another source of singularities in the distribution function in 3D. 
Indeed, j{9 1 r) becomes singular at 9 = or 9 = tt whenever sm(9) — 0, provided the corresponding 9q (called 9 g ) 
has a non-zero sinus-value. We call this singularity glory in analogy with the phenomenon well known in atmospheric 
optics and quantum-mechanical scattering theory [26,24]. As it will be shown latter, the classical singularities are 
replaced by pronounced peaks accompanied by fast oscillations in the quantum regime. 

As follows from Eq.(6), the glory appears at 9 = when the equation 

e g = Prsm(6 g ) 

has a non-trivial solution (see Fig. 3 (b)). The glory appears first exactly at the focusing time 17 = 1/P, and stays 
forever after that. The glory peak of the distribution at the north pole is shown at Fig. 4. Note that both 9q and 9 g 
are zero at Pt = 1 . The backward glory appears at 9 — tt after the rainbow ring reaches the south pole at t = r£ . For 
r > t'j , two glory angles 9' gl , 9' g2 exist for which 

n = 9' gs -Pr S m(9' gs ) s = 1,2 

After being created, the backward glory stays forever in the classical model. Figure 5 shows similar phenomena for 
the polarization-type interaction. 



C. Quantum kicked rotor (2D) 



We proceed now to the quantum description of the kicked rotor, starting with the simplest two-dimensional case. 
The dynamics of a quantum 2D rotor after its interaction with a strong laser pulse was described in detail in [19,37]. 
We reproduce briefly the relevant results here. 

The wave function of a 2D rotor can be expanded in terms of cigenfunctions of a free rotator 



^ +00 

*(0,t) = -= Cn(T)exp(in9) (10) 



+ OC 

2tt in 

(here we use the dimensionless units). In the absence of the field, the wave function is 

1 +°° 2 
*(0, T ) = -= V c„(O)expH^-r + m0) (11) 

V Z7T — ^ 

m— — 00 
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Even in this simplest case, the wave function (11) shows an extremely rich space-time dynamics. In particular, r) 
demonstrates a periodic behavior in time (quantum revivals) with the "quantum" period T rev = An (full revival), 
and a number of fractional revivals at t — p/s T rev (p and s are mutually prime numbers) [22]. 

The result of an extremely short pulse action on the quantum system can be obtained by neglecting the effect of the 
kinetic energy term in the Hamiltonian (1) during the pulse. The reduced Schrodingcr equation can be immediately 
solved, providing the following rule connecting the wave function before and after the kick: 

*(0,r fe + O) =cxp[-^y dr V(6,t)} *(0,T fe -O) (12) 

Here Tfc is the moment when the kick is applied, and integration is performed over the short kick duration. In the case 
of the dipole- type interaction, the exponent in Eq.(12) takes the form of exp[iPcos(6>)]. Using Eq.(lO), we derive the 
transformation rule for the coefficients c„: 

+ OC 

ra— — oo 

Here J n (P) is Bessel function of the nth order. If, for example, we apply a kick to a rotor being in its ground state 
at r = (i.e. ^{0, 0) = l/\/27r), then the following rotational wave packet is created: 

*W> T ) = E ^Jn{P)^M-i\r + in6) (13) 

ra— — oo 

Figure 6 shows time evolution of the probability density \^f(6, r)| 2 numerically calculated according to Eq.(13) for a 
relatively large kick strength (P = 85). We may identify there some distinct phenomena that were already discussed 
in the classical context above. First of all, the wave function shows an extreme narrowing in the region S«0 (Fig. 6 
(b) ). Clearly, this effect may be interpreted as a focusing at r/ = 1/P. Immediately after the focusing, two sharp 
singularities are formed in the distribution, which are moving with time. They may be identified with the rainbow- 
type catastrophe described above. The quantum nature of rotational motion replaces the classical singularities by 
sharp maxima with the Airy- like shape typical to the rainbow phenomenon (Fig. 6 (c)). We will discuss this in more 
detail in the following sections presenting the semiclassical description of the rotational wave packet dynamics. 

It should be stressed, however, that the long-time regime of the quantum motion is radically different from the 
classical one. Thus, contrary to the classical limit in which the singular caustics exist forever, they gradually disappear 
in the quantum case because of the overall dephasing of the initial rotational wave packet. Figures 6 (d) - 6 (i) show 
several examples of fractional foci and rainbows in the angular distribution, which is a purely quantum effect. 

The results for a 2D rotor with polarization-type interaction are very similar to the ones given above, except the 
obvious symmetry differences. Therefore, we proceed to the analysis of the three-dimensional effects now. 



D. Quantum kicked rotor (3D). 



The treatment of the three-dimensional case is rather similar to the previous one, and we will briefly summarize the 
main differences. First of all, for a vertically polarized (along z-axis) external field, the spherical harmonics Y[ (£l) 
with = (y>, 9) replace the eigenf unctions exp(in9) in the rotational wave packet superposition. Next, the energy 
spectrum is 1(1 + l)/2 instead of n 2 /2 now. A wave function of a molecule with permanent dipole moment undergoes 
the following transformation as the result of a <5-pulse applied at r = 0: 



V>(O,0 + ) = exp(iP cos 0)^,0") = 



£iV47r(2/ + l)j z (P)y J °(n) 



1=0 



V(fi,o-) 



(14) 



Here, ji(P) is spherical Bessel function of the order I, and ip(Cl, ),ip(Cl, + ) are the wave functions just before and 
just after the interaction, respectively. If the rotor is initially in the ground state, ip(£l,0~) = Y" ° = 1/VAtt, then the 
wave function at time r after the kick becomes 



V(0,t) = ]Ti z V2TTij ; (P) exp 



-l(l + l)r 



(15) 



5 



Fig. 4 (a), bottom line, shows the probability distribution \%P(8,t)\ 2 for a strongly kicked (P = 75) rotor at the 
focusing time t = 1/P (a), and also at r — 2/P (b), and t = 4/P (c). The presence of a glory at 9 = and of 
rainbow with Airy-like oscillations is clearly seen in these pictures. An analytical semiclassical description of these 
features will be given in the following sections of the paper. 

If the external field couples with induced polarization of the molecule (see Eq.(5)), the wave function immediately 
after the kick can be written as 



V>(ft,0+) = cxp (iPcos 2 9) V>(O,0") 



v^exp $>y» 



(16) 



with 



(17) 



The coefficients cLl.i can be determined by a recurrence relation (see Appendix A). For a rotor initially in the ground 
state, the wave function after a kick can be expressed as 



jL Vc/exp [-U(21 + 1)t}Y°(0). 
4tt Y 



(18) 



(Note that only even spherical harmonics contribute to the wave function here). Figure 5 shows the probability 
distribution \ip(9,T)\ 2 according to Eq.(18) for r = 0.5/ P (a), r = 2/P (b), and r = 2.4/ P (c), compared to the 
classical distribution function. As in the case of molecules with a permanent dipole moment, we observe focusing (a), 
rainbow and glory formation ((b) and (c)). 



III. SEMICLASSICAL APPROXIMATIONS 

We have observed a close correspondence between the classical angular distribution function and exact quantum 
mechanical probability distribution in the case of strong kicks. This invites to a quasiclassical treatment of the time- 
dynamics of a strongly kicked rotor. The following part of the paper presents an asymptotic quasiclassical theory 
of the above critical phenomena, which is based on integral representation of the rotor wave function with the help 
of quantum mechanical propagator. For simplicity, the discussion will be focused on rotors (molecules) having a 
permanent dipole moment. 



A. Propagator for a 2D rigid rotor 

The propagator for a free quantum mechanical 2D rigid rotor is 



^ oo 

U(0 o ,0;0,t) = -== 



p— — oo 



. (9 - 9 Q + 2tt P ) 2 
2t 



(19) 



The sum over p corresponds to an infinite number of trajectories going from 9q to 9 by making various number of full 
rotations. The wave function of the rotor is therefore 



^(6,t) = / <W o U(0 o ,O;e,T)il>(0o,O). 
Jo 

Since ip(9o + 2pir, 0) = V>(#o, 0), the wave function can be written as 



\ / 2l1TT J-, 



d9 n cxp 



. (fl-flo) 2 
1 2t 



V>(#o,0). 



(20) 



(21) 



If interaction with the laser field is proportional to cos# (see Eq.(2)), the wave function after the kick is 



G 



V(0,t) 



2ttVJt 



d6 exp 9 )} 



(22) 



with 



$(0,6o) = Pcos0 o + 



(g - ^o) 2 
2r 



(23) 



B. A planar model for a 3D rigid rotor 

Classical considerations about the focusing process and the formation of glory and rainbow suggest that those 
phenomena arise from particles that have been initially in the region of small angles (9 « 1). Therefore, we expect 
to get a certain insight by considering an almost planar motion of a rotor in the polar region. For this, we consider 
a fiat 2D model in which the wave function is initially uniformly distributed inside a circle of area 4tt and is kicked 
radially towards the center of the circle by a cosine kick. 

According to this model, the wave function at time r after the kick can be written as 



1 



where I is defined as 



with the phase 



Aim 3 / 2 



exp 



1=1 d9 d / dipo exp [i<$>(9, 6 , ip )] , 
Jo Jo 

<P(9, Oo, to) = + \ [l P) *8 ^ cos to- 



(24) 



(25) 



Since we assume that only small values of 9 are important, we expanded cos6* up to fourth order. Integrating 
Eq.(24) over </? , we find 



iTV4ir 



exp 



9 2 

P+ Yr 



o 



'UK, II, -I,, (^y- ) exp 



p 

24' 



(26) 



where Jo is Bessel function of zero-th order. Figure 7 shows the probability distribution \ip(9,T)\ 2 for a fully three- 
dimensional rotor system (Eq.(15)) compared to the probability distribution according to Eq.(26) at different times. 
It can be seen, that for small r, the planar model describes the dynamics of the rigid rotor quite accurately. Focusing 
(Fig. 7(a)) and glory formation (Figs. 7(b) and (c)) appear also in the planar model. Only at larger times (see Fig. 
7(d)), when the wave function of the 3D rotor depends on contributions from larger angles 9, the planar model fails. 
Moreover, by analytically calculating the wave function at the focal point (see Appendix B), we show that the wave 
function in the planar model has the right asymptotic behavior for large P. In Sec. IV. B, we will use the planar 
model exactly as it was defined above to describe the focusing of the wave packet which occurs at zero angle and for 
short time. In Sections V.B and VI, an extension of the model will be practiced since we will use the full P cos 9o 
term inside the phase. Although not justified for the large angles that we will consider, we feel that the physics of 
the planar model is sound enough to provide a qualitative understanding of the behavior of the wave packet on the 
sphere for not too large times. 



C. Topology of stationary points 

For a large kick strength P» 1, the exponent in Eq. (22) and Eq. (24) is a fast-oscillating function. Therefore, the 
main contribution to the integrals comes from the points of stationary phase, where d$/d9o = (and dQ/dtpo = 
in the 3D case). Figure 3 shows the points of stationary phase for every angle 9 for (a) 2D and (b) 3D rotors. The 
shape of the wave function is determined by the topology of the points of stationary phase. The simplest semiclassical 
approximation is the stationary phase approximation, which can be used if the stationary points are well separated. 
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However, the semiclassical catastrophes like rainbow and glory formation occur when the stationary points are close 
to each other and, even, merge. 

Thus, if Pt w 1, the three points of stationary phase around 9 w (and for 9 w 2ir in the 2D case) are non- isolated 
(cusp). For Pt — ► 1, they merge into a single one, giving rise to focusing. In this region, all three stationary points 
have to be treated together. In Sec. IV we show how the wave function can be approximated with the help of Pearcey 
approximation in this case. 

Next, we consider the case of Pt > 1. For 9 w r , we have two close stationary points merging into a single one 
at 9 = 9 r . This fold-type catastrophe causes the formation of angular rainbow which can be treated with the help 
of Airy approximation, as we will show in Sec. V. In three-dimensional case, the points of stationary phase belong 
to two different branches (see Fig. 3 (b)). The point with 9 > 9 g corresponds to ip = 0, while the two points with 
9 a < 9 g belong to the branch with <p = ir. These two branches merge at 9 = 9 g giving rise to the glory formation, 
which is absent in the 2D case. In Sec. VI we show how this phenomenon can be treated with the help of Bessel 
approximation. 



IV. FOCUSING 



A. Focusing of a 2D rotor 



In this section we consider 2D kicked rotor close to focusing ( Pt w 1). As can be seen in Fig. 3 (a), there are three 
close points of stationary phase near 9 w 0, and three stationary points at 9 w 2ir in this case. Due to the symmetry, 
it is sufficient to consider only the stationary points at 9 w 0. The complete wave function Eq.(22) can be written as 

1>(0,T)=j>{0,T)+$(2ir-6,T). (27) 
In order to account for the three stationary points, we keep all terms up to the order # 4 in the phase Eq.(23). Then, 
1 



V(0,t) 



: exp 



Making the following substitutions 



2t 



P 



— I 

s/2^ J-, 



exp < % 



24 2 



P 



d9 



(28) 



u = 



1/4 



#0, 



X = 



"6" 


1 

2 










P 







we can write the wave function Eq.(28) as 

4>(0,t) = 



1 


"6" 


1 
4 






TX\j2lT 


P 


exp 







where V is the Pearcey function 

/oo 
exp[i(u 4 + xu 2 + f3u)} 
-OO 



and (3 = V2 - 



V{x,(3) , 



du 



(29) 



(30) 



1 00 00 „m aln 

2 ^ ^ ml 2nl 

n—0 m—0 



-(2n + 2m+ 1) 

4 



exp 



i-(Wn + 6m + 1) 



This function has been studied intensively (see, e.g. [35,25]. If Pt = 1, that is, at focal time, then x = 0, and the 
wave function Eq.(29) can be reduced to a single sum, 



1 



277^ 

2n\ 



= (6P) 1 / 4 exp 



iP It 



00 aln 



n=0 



-(2n + l) 



exp 



i-(10n+l) 



(31) 



At the origin 6> = the probability distribution is 



8 



8tt 2 



<6PT Z - 



0.4078VP. 



(32) 



For P — > oo, the probability distribution at = goes to infinity, which is consistent with the classical picture. 
Moreover, we can determine from Eq.(31) that for large angles 9, the probability distribution is independent of P 



V> [0 



1 



7T(60) 2 /3- 



Considering the classical angular distribution function of a 2D rotor (see [19,37]), 

1 1 



f(0,l/P) = 



2tt \dO/de \ 



(33) 



(34) 



with 9{t) = 9 q -Pt sin(0g) (mod27r), we obtain the same result by expanding sin0g up to the cubic term. 

In Fig. 8 the wave function in the Pearcey approximation Eq.(30) is compared to the exact wave function. Fig. 8 
(a) demonstrates that the focusing is very accurately described in the Pearcey approximation. For r = 1.1 jP (b) and 
t = 1.2/ P (c), we still observe a good agreement with the exact quantum mechanical result, but for larger times (see 
t = 1.4/P in Fig. 8 (d) ), the approximation starts to break down. 



B. Focusing of a 3D rotor 



We proceed now to the case of Pt w 1 for a 3D rotor (see Fig. 3 (b)). In order to treat all three points of stationary 
phase around 9 — together, we have to keep the full quartic expression in the phase $(#, #o, <A)) m Eq.(25) to 
accurately calculate the integral. Due to the geometrical factor 9o in front of the exponent in integral Eq.(24), the 
integration over 9q can be expressed with help of a derivative of Pearcey integral. We can write the wave function as 



W,t) = - 



4 T7r 3/2 



exp 



i [P + 



9J_ 

'It 



1/2 ,.2ix 



dVi(x, y) 



dfo- ., 
o dy 



where 



/>OC 

^M^i y) = / rfwcxp [i (u 4 + xu 2 + yu)] 
Jo 



(35) 



(36) 



is a part of Pearcey integral, such that V(x, y) = Vi(x, y) + Vi(x, —y). Here, y = — /3cos<po- According to [36], the 
integral V\ can be represented as a sum: 



oo oo 



m n 

V\ = - — — r 

n—0 m—0 



i(n + 2m + l) 
4 



exp 



ii7r (n + 6m + 1) 



exp ( i— Tin 



Since 



L 



2 \ n-i o («-2)!! 
dtpo cos ip = 2tt- 



for odd n and zero if n is even, the wave function Eq.(35) can be written as 



(37) 



(38) 



6 



1/2 



p I CX P 



i,P+ 27 



1 ^ ^ x m f3 2n (2n- 1)!! 

n—0 m—0 



4^rf-f- m! (2n)! (2n)!! 



x r 



-(n + to + 1) 



exp 



7T "1 

i — (5n + 3to + 3)J . 



(39) 



Comparison of the results obtained with this expression with those resulting from the exact quantum evolution for 
Pt E [1,1-4] and 9 € [0,0.3] is presented in Fig. 9. As in the 2D case, we can see that Pearcey approximation 
describes the focusing and glory very well at times Pt w 1. All three stationary points merge into a single one 6>o = 
at focusing time Pt = 1. In this = 0, and Eq.(39) can be simplified to 
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(6P) 1/2 cxp 



* ,P+ 27 



E 



/3 2 ™ (2n-l)!! 



V^ (2n)! (2n)!! 



cxp 



i-(5n + 3) 



(40) 



Fig. 9 (a) shows the focused wave function according to Eq.(40) compared to the exact quantum result. At the origin 
= 0, wave function Eq.(40) is 



V o, 



p 



3P 



-P. 



(41) 



As in the 2D case, the probability distribution diverges as P — > oo. However, in the 3D case, the probability 
distribution is proportional to P , while in the 2D case, it is proportional to y/P. 

Comparing Fig. 8 (c) and (d) with Fig. 9 (c) and (d), one can recognize the difference between the 2D and 3D cases. 
In the former case the peak at — splits into a symmetric pair of peaks moving away from = with time. In the 
3D case, the peak stays at 9 = after being created. This is the result of the glory effect mentioned before in the 
classical context, and discussed in detail in one of the following sections. 



V. RAINBOW 



A. 2D rotational wave packet in terms of Airy functions 



In the following, we investigate the formation of angular rainbow, that is, the shape of wavefunction around the 
rainbow angle for Pt > 1. In the 2D case, we can approximate the wave packet Eq.(22) as a sum of two contributions 
related to the rainbow angles 9 — 9 r and 9 = 2tt — 9 r (see Fig. 3 (a)). Due to the symmetry we can write 

m r) « i>e r (0, t) + Ve r (2tt - 0, t) . (42) 

The rainbow angle 9 r is defined as 9 r — 9~ — Pt sin 0~ o , and the phase Eq.(23) can be expanded around 9 r as 

2+(0-0 o ) 2 



*(Mo) 



2t 



+ _ - 9o) + p sin g- W> -J o)3 + . . .. 



3! 



Using the definitions 



Psin O 



1/3 



(00 - 0o ) and r\ = 



Psin O 



1/3 



0) 



(43) 



(44) 



the phase <& can be approximated by 



$(0,0 O ) 



- 0o) 2 



2t 



Within this approximation 



1 



27r\/iT 



exp 



2 + (0 - 0"o) ; 
2r 



+ VV+ y- 



1/3 



Using formula Ai(r]) = ^ J °° cos(y 3 /3 + rjy) dy one obtains 



1 


2 


1/3 






Psin O _ 


exp 


H 


y/ir 







P sin O 



' 2 + (0-0" o ) ; 
2r 



»(w+y 3 /3) 



dy. 



Ai( V ). 



(45) 



Figure 10 illustrates the quality of approximating the 2D wave packet by a sum of two contributions expressed 
through Airy functions around rainbow angles Eq.(42). We have chosen Pt = 4,4.7 and 6 in (a), (b) and (c), 
respectively. The case (a) corresponds to the rainbow angle 9 r = 2.555 and the two wave structures moving on the 
circle in opposite directions are still well separated. The interference between them is so small that it is not seen in 
the figure. For case (b), r = 3.236 and the two semiclassical wave functions interfere strongly, reproducing quite well 
the structure of the exact quantum wave packet. In case (c), r = 4.513 and the interference region extends between 
7r/2 and 37r/2 (more precisely, between two rainbow angles). Again, the semiclassical wave function reproduces very 
well the qualitative features of the quantum wave packet. 
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B. Uniform Airy approximation for rainbow in 3D case 



In order to calculate the 3D wave function for Pt > 1 and 9 ss 9 r , we have to estimate the double integral Eq.(24). 
Using stationary phase approximation for the integration over ipo, we obtain 



-Z7TT 



/ d8 6> 
Jo 



exp 



(9 + 9 
2r 



P COS #0 



(46) 



Here we consider only the azimuthal stationary point ifo = ir that leads to two stationary points #03 and #02 in Fig. 3 
(b). The second stationary azimuthal point ipo = is related to the isolated stationary point O i (Fig. 3 (b)) and its 
contribution to the rainbow wave function is negligible. To treat this important 3D case, we rely on the uniform WKB 
approximation that is supposed to be superior compared to the plain asymptotic result based on expansion Eq.(43). 
The wave function in the uniform Airy approximation is given by 



1 



Aim 2 



I . 



(47) 



Since the phase Eq.(25) close to the stationary points 9 2 and #03 has almost cubic structure (see (43)), one can 
approximate I in the region 8 < 9 r as the integral 



J = 



-VKT 



29 



dx 



,9i + 92X 
2ir 



exp [i(A + £x + x 3 /3)] 



(48) 



We fix parameters <?i, (/2, A, and ^ in such a manner that the values of I and J coincide at the points where their 



phase is stationary. All these quantities depend on 8. The stationary points of I satisfy £ 
x\ = — and X2 = \/—%- From identification of the phases we finally obtain 



A = 



$(6,603) + $(0,M 



0. As C < 0, 



(49) 



and 



3 



($(fl,fl03)-*(^M) 



2/3 



(50) 



where 



*(0,0 Oa )= ( + e °* )2 +Pcose 0s , 

It 



(51) 



with s = 1,2. Identification of the integrals leads to two linear equations for gi and 52- Solving them, one obtains 



9i 



Hp] i«r 



?02 



I cos 9 — COS 9 2 I 



COS 9 — COS 003 I 



(52) 



and 



.92 = 7T 



0: 



02 



I COS 6*o — COS 6*02 I 



#03 



I cos — cos 6, 



03 I 



(53) 



For 9 = 9 r the coefficient (72 approaches zero. Hence, for > r the uniform Airy approximation reduces to Airy 
approximation 



*^) = iv^ Gie ^ Ai(??) 



(54) 



where 
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Gi = lim gi = 2tt 



Psin O 



(55) 



and rj is defined in (44). Calculations performed according to these formulas are presented in figure 11. The uniform 
Airy (UA) approximation gives quite good description of the quantum wave function. The figures do not show the 
wave function shape and value for very small where the UA approximation fails (singularity l/y/0). The quantum 
wave function is normalized (2tt f£ \^f(0, r)| 2 sin0d0 = 1). For the wave function in UA approximation this integral 
gives 0.838. The main difference between the Airy approximation for the 2D rotor and the UA approximation for the 
3D rotor is the l/y/0— factor in the 3D case, which leads to the divergence of the 3D semiclassical wave function at 
very small values of 9. This singularity, however, does not show up in the probability density as the latter is given by 
271-^(0, t)| 2 sin 0. 



VI. GLORY 



In the 3D case, for Pt ^ 1 and 9 ~ 0, we have an isolated stationary point around 9 = 0, and two close stationary 
points (#02 and #01) around the glory angle 9 g (see Fig. 3 (b)). The latter pair of stationary points is responsible for 
the formation of glory, and it can be treated by the uniform Bessel approximation [34] . According to [34] , we perform 
the integration over #0 in Eq.(24) by using the ordinary stationary phase approximation: 

(56) 



2V27 f #0 6o{tpo) ol/2 cxp [i me (<po)} - X)} 



|<9 2 $/50 2 | 

Here x = — tt/4 if <9 2 <I>/<90 2 > and \ = tt/A in the opposite case. Angle 0o(<£>o) stands for the stationary value of 0o 
for a given value of (po- In order to perform the remaining integration, we map the phase $ [0q(</?o)] with the function 



$(0,V) = a{9) - b{9) cosip. 
The mapping is one-to-one if <po = corresponds to ip = 0, and <£o = t to ip = ir. Then 

a{9) 



i[$(M + <lKM] and b{6) = X - [$(0„2) - • 



(57) 



(58) 



Here, 0oi and 002 are the stationary points corresponding to tp = and ip = ir, respectively. The integral in Eq.(56) 
can be written as 



I = 2V2^ : exp [i (a(0) - x)] / dtp 



/0 ' |(l_p r) + ^2|V 

By differentiating <J> with respect to ip, we obtain (using d&/d9 = 0), 

- — — = sin^o-77- = b{0) sm ip. 

oipo dip t dip 

By taking another derivative of Eq.(60) at ip — and ip — ir we arrive at 

difo 



Y^-exp [-ib(0) cosV'] ■ 



which determines d(po/dip at the stationary points. Following the derivation of Berry [34], we define 

= p + (0) + p^ (0) cos ip. 



Q dipo 



|(i-Pr) + ^0 o 2 | 1/2 

According to the mapping conditions, p± can be determined as 

\ V2 



(59) 



(60) 



(61) 



(62) 



P±(0) = 2' 



1/2- 



7{)1 



(l_p T ) + ^2 i 



± 



702 



( l_p r) + ^2 



(63) 
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The integral can then be written as 



I = 2^2^57 exp [i (a(9) - x)] [ # [p+(0) +P-(0) cost/;] exp [ib{6) cos-0] 

Jo 

- 2^\/2^cxp [i (a(6) - x )\ {p+(0)Jo [6(0)] - ip_(fl)Ji [&(*)]} • 



The wave function is 



(2tt) 3 / 2 t\/2 



exp 



J. 



(64) 



(65) 



Figure 12 shows the modulus squared of the wave function (in the uniform Bcssel approximation) for P = 75 and 
t = 4/P in comparison with the exact numerical integration of Eq.(26). For small angles 9, the quantum mechanical 
wave function is well represented by the uniform Besscl approximation, and the formation of the angular glory at 
9 w is reproduced. However, the approximation breaks down if 9 approaches the rainbow angle 9 r . 

If 9 — > 0, the stationary points O i an d $02 merge to a single stationary point 9 g . For time only slightly larger than 
the focusing time, 9 g is defined by 



g = 



Pt-1 
Pt 



(66) 



and the wave function becomes 
1 



W2t 



exp[i(a(0) - x)] 



Mo (¥) 



exp 



2$ 



(67) 



This expression is analogous to Ford- Wheeler formula in nuclear scattering [26]. It can be also derived by directly 
calculating the integral in Eq.(26) with the help of simple stationary phase approximation using the glory angle 9 g as 
a single stationary point. 



VII. KICKED 3D ROTOR AT FINITE TEMPERATURE (CLASSICAL TREATMENT) 

In the previous sections, we considered a quantum rotor being initially in its ground state, or, in the classical 
context, an ensemble of rotors with zero initial velocity. In the present section we investigate the dynamics of an 
ensemble of kicked classical rotors with a finite initial temperature, and compare the results with the zero-temperature 
case. 

The time evolution of a 2D rotor under thermal conditions was studied in [37]. In the previous Sections we 
have already observed several differences between 2D and 3D kicked rotors at zero initial temperature. Another 
difference appears at finite temperature because of the conservation of the angular momentum projection onto the 
field polarization direction. In a thermal ensemble, molecules generally have a non-zero initial projection of this kind. 
As a result, an effective repulsive centrifugal force prevents the rotors from reaching the exact 9 = and 9 = tt 
orientations. Therefore, two holes in the angular distribution of a driven 3D thermal ensemble should be always 
present at 9 = and 9 = n (sec a related discussion in [41]). 

In the following, we describe the dynamics of a classical ensemble of 3D rotors under thermal conditions. The 
equation of motion for the angle 9 of a free 3D rotor is given by 

m = §^, (68) 

l z sin 9 

where the canonical momentum p v is a constant of motion. We introduce dimcnsionlcss variables by defining the 
thermal momentum p t h = I^th with uj t h = y/ksT/I, where T is the temperature and ks the Boltzman constant. 
Equation (68) can then be written as 

sin 9 
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where p' = p v /pth is the dimensionless momentum and t' = Lo t ht is the dimensionless time. The differential equation 
Eq.(69) can be integrated by substituting u = cos8, so that 



cos 



0(0 = l (l - — ) cos [6(0) + tot'] + + cos [6(0) - tot'] . 

2 y LU J 2 y LO J 



Here. 



1/2 



■LJ 



sin 2 6(0) 



If a kick is applied to the rotor at t' = 0, it gains the momentum 

p' e =p' e (0)-P'sm6(0). 



(70) 



(71) 



(72) 



Here, p' g (0) = Pe(0)/p t h is the initial canonical momentum in dimensionless units, and P' — ^l[ j dt£(t)\/ (Iuj t h) is the 
dimensionless kick-strength. If P' m 1, the energy gained due to the kick is comparable to the initial thermal energy 
of the rotor. 

At zero initial temperature, p'^0) — p' e (0) — 0. In this case, Eq.(70) is reduced to 



0(0 + 0(0) -ft 7 sin 6(0). 



(73) 



Note, that while P' — > oo for zero temperature, the product P't' remains finite and corresponds to the product Pt in 
Eq.(6). 

At finite temperature T, the initial distribution function of the rotors is 



f{Q,<P,p'e,P' v ) = q exp 



(74) 



with Q being the partition function. Note, the p' g and p' v depend on the temperature. 

In the following, we present the results of Monte-Carlo simulation of the dynamics of a thermal ensemble of kicked 
rotors. Figure 13 shows the time evolution of the angular distribution of an ensemble of molecules with permanent 
dipoles subject to a kick of strength P' = 10 (solid line), P' = 5 (dashed line), and P' = 1 (dotted line). In (a), the 
angular distribution is plotted at focusing time P't' = 1. For a strong kick, we can still observe focusing, although 
the singularities in the zero-temperature case are replaced by pronounced peaks. Moreover, the "repulsive centrifugal 
force" causes a hole around = 0. The widths of the hole is increasing with increasing the temperature (or decreasing 
P'). The angular distribution for P't' = 3 and P't' — 4.5 is shown in Fig. 13 (b) and (c), respectively. For P' = 10 
and P 1 = 5 the remnants of the angular rainbow moving towards the south pole are still seen. If the initial thermal 
momentum of the particles is comparable with the momentum gained through the kick (P' = 1), the kick results only 
in a slight perturbation of the equilibrium distribution (see Fig. 13 (a) - (c)). 



VIII. MULTIPLE KICKS: ACCUMULATIVE SQUEEZING OF THE ROTOR ANGULAR 

DISTRIBUTION. 

In the previous chapters, we investigated the wave function of a rigid rotor kicked by a single ultra-short laser pulse. 
We saw that the initially uniform angular distribution becomes much narrower in angle after the kick. However, the 
degree of orientation (or alignment) that can be achieved by a single laser kick is limited. It has been shown by us 
recently that the angular distribution of an ensemble of rotors can be further squeezed by applying a series of short 
laser pulses. A generic approach, the accumulative squeezing strategy, was proposed in [19]. According to this strategy, 
the time-dependent angular distribution of the rotors is measured after each pulse, and the next pulse is applied at 
the time of the minimal angular spread. It was shown that for 2D rotors the width of the angular distribution is 
gradually decreasing with number of kicks without any sign of saturation [37]. In principle, an "unlimited" angular 
squeezing may be achieved this way. 

In the following, we show how accumulative squeezing works for 3D rotors - with permanent dipole as well as for 
interaction with induced polarization. We also study analytically the asymptotic behavior of accumulative squeezing 
after a large number of kicks and show, that the squeezing is indeed unlimited in this approach. As a measure 
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of orientation we choose the orientation factor defined as O(r) = (1 — cos#(t)). in a similar way, alignment is 
characterized by the alignment factor A(t) = (1 — cos 2 0(t)). 

We first provide results for angular squeezing of an ensemble of classical rotors. Figure 14 shows the orientation 
factor (alignment factor) as the function of the number of applied short pulses. The dots represent an ensemble 
with initially zero temperature, and the crosses represent two ensembles with different finite initial temperature. In 
all three cases we observe that angular squeezing is gradually improving by applying more and more kicks, without 
saturation. 

As the next step, we investigate analytically the squeezing process after a large number of kicks, i.e. when the 
angular distribution is already very narrow. The most of the particles are close to the pole region at this stage, so 
that we can neglect the spherical geometry and consider a flat model. Moreover, since the distribution is uniform in 
the ^—coordinate, we can transform the problem to dimensionless Cartesian coordinates x — 8 cos ip and y — 9 sirup, 
and consider a distribution having the same width in x- and y— directions. Because the interaction with the field is 
almost harmonic in the polar region, we can separate the coordinates and consider independently the behavior of the 
system spread in x and y coordinates. 

Therefore, let us focus on the one-dimensional accumulative squeezing of a particle moving along x-axis. The 
classical and quantum analysis of this situation are very similar to each other, so we provide the quantum version 
here. We assume that 



x(t) = p(r) = 0. 



(75) 



because of the symmetry. For freely moving particles, the spatial dispersion x 2 (t + At) at time r + At depends on 
three distribution moments x 2 (t), P 2 (t), and xp + px(r) defined at the initial time r: 



x 2 (t + At) = x 2 {t) + xp + px(t)At + p 2 Ar 2 . 



In the similar way: 



xp + px(r + At) = xp + px(r) + 2p 2 (t)At. 



(76) 



(77) 



It follows immediately from Eq.(76) that r k is the time moment of extreme squeezing if the mixed moment xp + px(r k ) 
is zero. Assume now that the system is kicked exactly at that instant of time (i.e., the wave function is multiplied by 
a factor exp (— iPx 2 /2)). The x 2 (r k ) remains unchanged, however the mixed moment takes a new value 

xp + px{T k ) = -2Px^{ Tk ). (78) 

Since it is negative, the modified wave packet will reach a new minimum in width at r k+ i = r k + Ar k , where 

Px 1 ^) 



p 2 (r k ) + Px 2 (r k ) 



(see Eq.(77)). There is a gain in squeezing 

x 2 (T k+1 ) = x 2 (r k ) 



1 - =- 



Px 2 (r k ) 



p 2 (T k ) + Px 2 (r k ) 



but the kick has increased the spread of momentum since 

p 2 -(T k+1 )=p 2 -(T k )+Px 2 ~(T k ). 

We introduce the new notations 

x 2 (r k ) = u k 
pP 2 (n) = w k . 

The quantities u and w satisfy the finite-difference equation 



(79) 



(80) 



(81) 



(82) 



u k+ i =u k ■ 

u k + w k 

W k+ l = w k + u k . 



(83) 
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For large k they can be replaced by differential equations 

du u 2 



dk w + u 

dk= U - (84) 

We notice that u = x 2 is monotonically decreasing with k, while w, that is proportional to the average kinetic energy 
is gradually increasing. The above system of differential equations may be solved exactly as it has an integral of 
motion 

u + 2wu = const. (85) 

However, in the regime of strongly developed squeezing (u — > 0), the differential equations may be further simplified 
to 

du u 2 

dk w 
dw 

dk = U ™ 

Then, the integral of motion is 

wu = const., (87) 
which is, of course, consistent with Eq.(85) in the limit u <C w. We finally arrive to the equation 

S=- c « 3 - < 88 » 

from which one immediately obtains the asymptotic behavior 

x 2 {r k ) ~ u ~ 



W 2 (Tk) ~ — ~ Vk 
u 

u 1 

w k 



u 

V 1 

At* ~ - ~ -. (89) 



IX. SUMMARY 



In this paper we concentrated on the semiclassical time evolution of two-dimensional and three-dimensional kicked 
rotational systems. In the 2D case we have developed a semiclassical description of the dynamics starting from the 
exact propagator. In the 3D case, the problem can be greatly simplified by approximate mapping the semiclassical 
motion on a sphere to the motion on a plane. Since we are studying wave packets in the vicinity of the poles of a 
sphere, this approximation is rather justified. 

We have shown that the quantum mechanical evolution of kicked rotors can be understood in the frame of the theory 
of catastrophes, a well known subject of mathematical physics. The kicked rotational systems experience either two 
or three catastrophes with time, depending on the dimension. Similar catastrophes have been studied in different 
context in various physical systems in the past. An interesting property of the kicked rotor is that the catastrophes 
appear in the time domain and can be sharpened by simply increasing the strength of the kick. 

The formation of a cusp is the first dramatic event in the dynamics of a kicked rotor, leading to the formation of a 
sharply focused rotational wave packet. In the classical limit, the cusp in the angular distribution exists for arbitrary 
small kick strength. In the two-dimensional case, the focused wave packet is very accurately described in terms of 
Pearcey's function. Due to the additional azimuthal coordinate, the wave packet in the three-dimensional case is 
expressed in terms of a modified Pearcey function. 

Just after the cusp, a fold singularity is created in the wave function, similar to the rainbow formation in scattering 
theory. The shape of the wave packet in the rainbow domain is explained in terms of Airy's function or with the 
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uniform Airy approximation. For two-dimensional rotors, there are two rainbows in the angular distribution which 
travel in the opposite directions and interfere with each other. In the three-dimensional case, another spectacular 
catastrophe occurs besides the rainbow: the glory effect. Semiclassically it is explained in terms of classical trajectories 
coalescing at 6 = and n, and the shape of the singular wave packet can be well described by the uniform Bessel 
approximation. 

In addition to zero-temperature results, we also considered the dynamics of kicked thermal ensembles of classical 
three-dimensional rotors. We demonstrated that thermal effects lead to holes in the focal and glory peaks in the 
angular distribution. Finally we have investigated analytically the asymptotic behavior of the angular squeezing 
process in the limit of large number of kicks. We proved that there is no saturation in this process. 

Our theory for squeezing of two-dimensional kicked rotors has already been used to optimize the squeezing of atoms 
in a pulsed optical lattice [36,37]. The orientation or alignment of molecules with electromagnetic pulses is another 
subject which can benefit from our study, the optimal conditions for molecular alignment by trains of short laser 
pulses have been recently studied in [14]. 
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APPENDIX A: RECURRENCE 



In this Appendix, we want to express exp (iP cos 2 9) in terms of spherical harmonics. We start with writing 

P P 



1 ' 1 + ~2 COs29 



exp (iP cos 2 6) = cxp 

= exp (if)E^ (f ) (2/ + 1) (cos 20) 

We then expand Pi (cos 29) in terms of Pi (cos#), 

P ; (cos20) -^d L , z P L (cos0). 



(Al) 



(A2) 



We apply the recurrence formula for Legendre polynomials to Pi (cos 29) and then expand Pi (cos 29) in terms of 
Pi (cos#) with cos 29 = u = 2x 2 — 1 and cos# = x. 



(I + l)Pi+i(u) = (21 + l)uPi(u) - Wi-i{u). 
(I + 1) d L .i+iP L (x) = (21 + l)(2x 2 - 1) 2 <fe,J^(s) - 1 Yl d L ,i-iP L (x). 



(A3) 

L L L 

Then, we multiply this equation with Pl> , integrate over x and use the orthogonality relation for Legendre polynomials 

2 



P L (x)P L ,(x)dx 



2L+1 



?L.L> 



and arrive at 



(l + l)d L j +1 ^-j = (2l + l) (^2 d ",i J 2x 2 P L (x)P L ,{x)dx 



2L + 1 



dh.i 
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2L+1 



dh,. 



(A4) 



(A5) 



Then, we introduce the (symmetric) matrix Nl^i = J Pl(x)Pl> (x)x 2 dx which has the matrix elements 



Nr. J, = 



(2L + 1) 2 
N L . L , = 



'{L + lf L 2 

2L + 3 + 2L-1 
2(L + l)(L + 2) 



N r . r, = 



(2L + l)(2L + 3)(2£ + 5) 
2L(L- 1) 



for 


L' 


= L 


for 


L' 


= L + 2 


for 


L' 


= L-2 



(2i + l)(2L + l) 
d L , +1 - — 



21 + 1 j I, 



i " L ' (2L-3)(2L-1)(2L + 1) 
and A^j,/ = otherwise. The recurrence relation is then 

^d L >.iN L . L , 

. L' 

From the properties of Nl,l> follows the selection rule dh,i = if L is odd. Furthermore, since Pl(1) = 1, 

L 

Note, however, that 7^ if Z is odd. 



(A6) 



(A7) 



(A8) 
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APPENDIX B: ASYMPTOTICS FOR THE WAVE FUNCTION IN THE FOCAL POINT FOR THE 

PLANAR MODEL 



Considering an ensemble of classical rotors, we know that at focal time all particles for which sin^o ~ are at 
the (north-) pole of the sphere. We want to investigate the semiclassical wave function in order to see if these 
considerations are also valid in the semiclassical approximation. Therefore, we calculate the wave function at focal 
time Pt = 1 at 6 = 0: 



is/Air 



cxp 



with 



I = P 



[ d8 8 
Jo 



cxp 



i— 61 
24 



(Bl) 



(B2) 



where L is the radius of the circle. Expressing the exponential function as a power series, we can perform the 
integration and obtain 



2 ^ v\ 2v+l 

v=0 



With 



(2i/)! 



2u + l 2, + !)! (I); 
we can express the sum as a confluent Hypergeometric function, 



PL 2 (1 3 .P 4 
/= — f 1 2'2' Z 24 L 



(B3) 



(B4) 



(B5) 



Now, we can investigate the asymptotic behavior. According to [43], for |PL 4 /24| 2 — > oo, the integral is approximately 



J w -i 



P 



m+ mp c ^\ l 2i L 



The absolute square of the wave function is therefore 

2 



V o 



i 



3P 
8tt 



24 



7T 



+ 



24 



PL 4 \PL 4 



1/2 



COS 



p 

24 J 



(B6) 



(B7) 



We can see, that for P — > oo, the function goes to the constant value 3P/8, which is independent of the radius L. 
This is in accordance to the classical results, where we find a singularity at = 0. 
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FIGURE CAPTIONS 

FIG. 1. 2D kicked rotor. Classical map representing the final angle 6 as a function of initial angle 9o for (a) r = 0.5r/, 
(b)r = Tf, (c) t = 3r/, and (d) r = 10t/. 

FIG. 2. Angular distribution of the initially uniform classical ensemble of 2D rotors after a single 5-kick. Figure (a) shows 
focusing at t — Tf, and (b) displays rainbows at r = 1.84r/. 

FIG. 3. Classical maps and critical angles for (a) 2D, and (b) 3D rotors. In the 2D case, the final angle is plotted modulo 
2ir in order to fold it into the interval [0, 2n]. In (a) Pt = 3, and in (b) Pt = 1.4. 

FIG. 4. Angular distribution of a linear molecule (rigid rotor) with permanent dipole moment. In (a), (b), and (c) the 
probability distribution is shown for Pt — 1, Pt — 2, and Pt = 4, respectively. In the top line, the probability of finding a 
classical rotor inside the solid angle element sin 9d9dip is plotted on a sphere. The middle line shows the classical distribution 
function as a function of 9, and the bottom line displays the corresponding quantum mechanical probability distribution for a 
kick-strength of P = 75. 

FIG. 5. Angular distribution of a linear molecule (rigid rotor) interacting with the external field due to induced polarization. 
In (a), (b), and (c) the probability distribution is shown for Pt = 0.5, Pt = 2, and Pt = 2.4, respectively. In the top line, 
the probability of finding a classical rotor inside the solid angle element sm 9 d9 dtp is plotted on a sphere. The middle line 
shows the classical distribution function as a function of 9, and the bottom line displays the corresponding quantum mechanical 
probability distribution for a kick-strength of P = 75. 

FIG. 6. Angular distribution of a 2D quantum rotor excited by a strong 5-kick (P = 85). The graphs correspond to (a) 
r = 0.5t/, (b) t = T f , (c) t = 2r/, (d) r = tj + T rev /2, (e) r = T f + T rev /3, (f) r = tj + T rev /4, (g) r = 1.8t/ + T rev /2, (h) 
r = 1.8r/ + T r ev/3, and (i) r = 1.8r/ + T re „/4, respectively. 



FIG. 7. Probability distribution in 9 for the 3D kicked rotor (solid line) and probability distribution for the planar model 
(dashed line). The figure shows the probability distribution for r = 1/P (a), r = 1.2/ P (b), r = 2/P (c), and r = 4/P (d). 
The kick strength is P — 50. 

FIG. 8. Figure (a) shows the probability distribution (2D case) at focal time r = 1/P. In (b), (c), and (d), r = 1.1/ P, 
t — 1.2/P, and r = 1A/P , respectively. The kick strength is P = 50. The solid line represents the exact quantum mechanical 
probability distribution, while the dashed line described the probability distribution in Pearcey approximation. 

FIG. 9. Glory and oscillations close to the focusing in the kicked 3D rotor. Figure (a) shows the probability distribution at 
focal time r = 1/P. In (b), (c), and (d), r = 1.1/P, r = 1.2/P, and r = 1.4/P, respectively. The kick strength is P = 50. 
The solid line represents the exact quantum mechanical probability distribution, while the dashed line shows the probability 
distribution obtained in the approximation given by Eq. (39) based on Pearcey function. 

FIG. 10. Airy approximation for 2D WP evolution for times (a) rP=4, (b) 4.7 and (c) 6, P—75. The solid line represents 
exact quantum result, the dashed one shows the Airy approximation. The arrow point to the rainbow angle 9 r . 

FIG. 11. Uniform Airy approximation in the 3D case (Eq. (47)) is compared to quantum result for Pr=4 and P=75. The 
arrow points to the rainbow angle. Figure (a) shows the probability density as a function of 9, and figures (b) and (c) show 
the real and the imaginary part of the wavefunction. 

FIG. 12. Absolute square of the wave function in uniform Bessel approximation (solid line) compared to the exact wave 
function (dashed line) for Pt = 4 and P = 75 (3D case). 
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FIG. 13. Angular distribution of an ensemble of classical rotors (in 3D) with finite initial temperature for P't' = 1 (a), 
P't' = 3 (b), and P't' = 4.5 (c). The solid line corresponds to the kick strength P' = 10, the dashed line to P' = 5 and the 
dotted line to P' = 1. 

FIG. 14. Accumulative squeezing of a classical ensemble of 3D rotors. Figure (a) shows the minimal orientation factor, and 
figure (b) shows the minimal alignment factor as a function of the number of kicks. The dots (.) correspond to an ensemble of 
rotors with zero initial temperature (P' — oo). The crosses (X) and (+) correspond to P' = 5 and P' = 1, respectively. The 
orientation factor and alignment factor were calculated with help of Monte Carlo simulation with 10000 particles 
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This figure "sc_Fig4.jpg" is available in "jpg" format from: 



http ://arXiv . org/ps/quant-ph/03 06 1 46 v 1 



This figure "sc_Fig5.jpg" is available in "jpg" format from: 



http ://arXiv . org/ps/quant-ph/03 06 1 46 v 1 




angle 



FIG. 2. Angular distribution of the initially uniform classical ensemble of 2D rotors after a single 5-kick. Figure (a) shows 
focusing at r = 77, and (b) displays rainbows at r = 1.84r/. 
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FIG. 3. Classical maps and critical angles for (a) 2D, and (b) 3D rotors. In the 2D case, the final angle is plotted modulo 
2-k in order to fold it into the interval [0, 2w]. In (a) Pt = 3, and in (b) Pt = 1.4. 
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FIG. 8. Figure (a) shows the probability distribution (2D case) at focal time r = 1/-P. In (b), (c), and (d), t = 1.1/P, 
t = 1.2/P, and r = 1.4/F, respectively. The kick strength is P = 50. The solid line represents the exact quantum mechanical 
probability distribution, while the dashed line described the probability distribution in Pearcey approximation. 
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FIG. 9. Glory and oscillations close to the focusing in the kicked 3D rotor. Figure (a) shows the probability distribution at 
focal time r = 1/P. In (b), (c), and (d), r = 1.1/P, r = 1.2/F, and r = 1.4/P, respectively. The kick strength is P = 50. 
The solid line represents the exact quantum mechanical probability distribution, while the dashed line shows the probability 
distribution obtained in the approximation given by Eq. (39) based on Pearcey function. 
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FIG. 11. Uniform Airy approximation in the 3D case (Eq. (47)) is compared to quantum result for Pr=4 and P=75. The 
arrow points to the rainbow angle. Figure (a) shows the probability density as a function of 6, and figures (b) and (c) show 
the real and the imaginary part of the wavefunction. 
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FIG. 12. Absolute square of the wave function in uniform Bessel approximation (solid line) compared to the exact wave 
function (dashed line) for Pt = 4 and P = 75 (3D case). 
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FIG. 13. Angular distribution of an ensemble of classical rotors (in 3D) with finite initial temperature for P't' = 1 (a), 
P't' = 3 (b), and P't' = 4.5 (c). The solid line corresponds to the kick strength P' = 10, the dashed line to P' = 5 and the 
dotted line to P' = 1. 



33 




number of kicks 

FIG. 14. Accumulative squeezing of a classical ensemble of 3D rotors. Figure (a) shows the minimal orientation factor, and 
figure (b) shows the minimal alignment factor as a function of the number of kicks. The dots (.) correspond to an ensemble of 
rotors with zero initial temperature (P' = oo). The crosses (X) and (+) correspond to P' = 5 and P' = 1, respectively. The 
orientation factor and alignment factor were calculated with help of Monte Carlo simulation with 10000 particles 
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